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ABSTRACT 

We present the spectral atlas of sources observed in low resolution with the Infrared Spectro- 
graph on board the Spitzer Space Telescope. More than 11000 distinct sources were extracted 
using a dedicated algorithm based on the SMART software with an optimal extraction (AdOpt 
package). These correspond to all 13 000 low resolution observations of fixed objects (both single 
source and cluster observations). The pipeline includes image cleaning, individual exposure com- 
bination, and background subtraction. A particular attention is given to bad pixel and outlier 
rejection at the image and spectra levels. Most sources are spatially unresolved so that optimal 
extraction reaches the highest possible signal-to-noise ratio. For all sources, an alternative ex- 
traction is also provided that accounts for all of the source flux within the aperture. CASSIS 
provides publishable quality spectra through an online database together with several important 
diagnostics, such as the source spatial extent and a quantitative measure of detection level. An- 
cillary data such as available spectroscopic redshifts are also provided. The database interface 
will eventually provide various ways to interact with the spectra, such as on-the-fly measurements 
of spectral features or comparisons among spectra. 

Subject headings: Atlases, Catalogs, Infrared: general. Methods: data analysis. Techniques: spectroscopic 



1. Introduction 

The Infrared Spectrograph (IRS; Houck et al. 
2004ill is one of three instruments on board the 
Spitzer Space Telescope (Werner et al. 2004) 
along with the two photometers Infrared Array 
Camera (IRAC; Fazio et al. 2004) and Multiband 
Imaging Photometer for Spitzer (MIPS; Rieke et 
al. 2004). The IRS performed more than 21000 
observations over the cryogenic mission lifetime 
(2003 November 30 - 2009 May 15), correspond- 
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ing to about 14000 distinct targets. The IRS 
observed between Ri 5.2 and w 38.0 /zm in two 
low-resolution modules with a resolving power of 
R ^ 60 — 120 (w 75% of the observations), and two 
high-resolution modules with R ^ 600 (Tabled]) . 

Two products are essential to ensure the legacy 
of the Spitzer/IRS data. Most importantly, the 
community should have available the highest qual- 
ity spectra for the simplest cases, i.e., isolated 
point-like or partially extended sources with neg- 
ligible (or uniform) background emission. In addi- 
tion, specific tools are required to analyze the most 
complex cases, i.e, extremely faint sources (few 
tenths of mJy), blended sources, and/or sources 
with non-uniform background emission. 

While existing data reduction and analysis 
software such as SMART (Spectroscopic Mod- 
eling Analysis and Reduction Tool; Higdon et 
al. 2004) and SPICE (Spitzer IRS Custom Ex- 
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tractioio) already provide several extraction algo- 
rithms, the AdOpt algorithm which we developed 
within SMART significantly improves spectral ex- 
tractions with the IRS by using a super-sampled 
point spread function (PSF; Lebouteiller et al. 
2010). AdOpt provides the best possible signal- 
to- noise ratio (S/N) while being a flexible algo- 
rithm able to handle blended sources and complex 
background emission. One of the first applications 
was the separation of a supernova and a galaxy 
nucleus located in the same slit (Fox et al (2010). 

The Cornell Atlas of Spitzer/IRS Sources 
(CASSIS) utilizes an automatic spectral extrac- 
tion tool based on AdOpt to provide the com- 
plete sample of staring observations (not including 
mappings) with the two low-resolution modules 
Short-Low and Long-Low (hereafter SL and LL). 
CASSIS complements the existing "post-BCD" 
databascll from the Spitzer Science Center (SSC), 
by providing optimal extractions as well as mul- 
tiple diagnostic tools including source extent and 
extractions with alternative backgrounds. Special 
attention is given to the identification/removal of 
artifacts, and to accurate background subtraction. 
The latter is especially important for faint sources, 
whose flux densities represent only a few percent 
of the total emission (mostly originating from zo- 
diacal dust at IRS wavelengths) . In the future the 
CASSIS atlas will include serendipitous sources 
as well as high-resolution spectra taken with the 
Short-High and Long-High modules. 

A dedicated website has been developecQ to 
allow users to query for observations and search 
for spectra described by various parameters. The 
CASSIS atlas currently contains about 13 000 low- 
resolution spectra corresponding to '--^ 11 000 dis- 
tinct sources. Table [2] shows the number of ob- 
servations in the CASSIS atlas for each scientific 
category. The first application of the CASSIS at- 
las is the study of 301 galaxies observed with the 
IRS and with the Infrared Astronomical Satellite 
(IRAS) (Sarsgyanet al. 2011). The CASSIS spec- 
tra were used to measure the 7.7 fim PAH flux den- 
sity and the warm dust continuum to derive the 
infrared luminosity. 



In the present paper, we review the most impor- 
tant steps of the extraction pipeline which leads to 
CASSIS (Sections 2 to 6). Complete and updated 
documentation is accessible onlin^. Products, di- 
agnostics, and the online interface are described 
in Sections 7, 8, and 9. Finally, we present a few 
examples and illustrations of the possibilities en- 
abled by the CASSIS atlas (Section 10). 

2. General description of the data process- 
ing 

About 85% of the IRS observations were per- 
formed in the "staring" mode which observes the 
spectrum of a given source at two different nod po- 
sitions along the sliiH. The two nod positions are 
sufficiently far apart that subtracting one nod im- 
age from the other removes the background effec- 
tively, except in cases when a serendipitous source 
is located at the other nod position, or when there 
is complex background emission. Note that back- 
grounds can also be subtracted using detector im- 
ages when the source is in the slit for the al- 
ternative order (e.g., nods in the LLl slit when 
the source is in the LL2 slit produce LLl images 
of background only). More details on the avail- 
able background subtraction methods are given 
in Sect. 13.3] In staring mode, the observation se- 
quence is performed as follows: 
order 2 nod 1 - order 2 nod 2 - order 1 nod 1 - 
order 1 nod 2 

For a given observation (called AORkey, or 
AOR; see the list of acronyms in the Appendix), 
the staring mode observes either a single position 
or multiple positions (the latter case being also 
referred to as "cluster observation" which is desig- 
nated by a single AORkey). The remaining < 15% 
observations were performed as spectral mappings 
in which the slit is progressively stepped across a 
extended source. We do not as yet include map- 
pings in CASSIS. 

The official SSC pipeline consists of two parts, 
the Basic Calibrated Data (BCD) pipeline and the 
post-BCD pipeline. The BCD pipeline reduces 
the raw detector images from the individual ex- 
posures and removes the electronic and optical ar- 



^http: / /irsa.ipac.caltech.edu/data/SPITZER/docs/dataanalysistools /toolsy^gpicey^ 
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tifacts, including dark current, droop effect, non- 
linearity, radhit (detection and flagging of cosmic 
ray events), jail-bar pattern, and stray light. The 
BCD pipeline also includes a flat-field correction. 
The steps performed by the BCD pipeline are rel- 
atively well understood, and we consider the BCD 
products as reliable inputs for spectral extraction. 
We refer to the IRS Instrument Handbook for 
more details on the BCD products. 

The post-BCD processing, which transforms 
calibrated 2D detector images to ID spectra, of- 
ten requires a manual intervention or an algorithm 
specifically designed for a given observation. The 
steps involved in the post-BCD processing notably 
includes image cleaning, background subtraction, 
and spectral extraction. 

The CASSIS pipeline provides an alternative 
to the post-BCD pipeline by enabling a flexible, 
automatic extraction algorithm that is able to 
handle very different observations, from barely 
detected sources to bright sources, from point-like 
sources to extended sources. CASSIS was built 
around the SMART- AdOpt tool (Lebouteiller et 
al. 2010) and uses several of its diagnostic tools 
(e.g., source extent, detection level). The basic 
steps of the pipeline are shown in Figures[l] and[2] 
and are reviewed briefly in the next sections; they 
can be summarized as follows: 

- Steps performed on the calibrated 2D im- 
age product: cleaning of bad pixels, co-addition 
of individual exposures, removal of background 
emission, 

- Spectral extraction: for a given source, two 
spectra are extracted, one for each nod position. 

- Steps performed on the ID spectra: com- 
bination of the nod spectra, flux calibration, de- 
fringing. 

The spectra are gathered in an online database 
which provides the single best overall spectrum, 
combined from both nods, with the maximum 
S/N. The database can also be queried to under- 
stand in more detail the various important diag- 
nostic parameters from the pipeline such as the 
inferred source extent or the detection level for al- 
ternative spectra. 



3. Image pre-processing 

The starting points of the CASSIS pipeline are 
the BCD images, more exactly the single expo- 
sure images also referred to as DCE (Data Collec- 
tion Events) images. Exposure images are carried 
through the extraction pipeline along with the cor- 
responding uncertainty image file and the mask of 
flagged pixels (BMASK, see Sect.[3D). 

The first public version of CASSIS uses data 
created with the SSC pipeline release S18.7.0. In- 
ternally, this version corresponds to the 4**^ iter- 
ation of the CASSIS pipeline, so that products 
are labelled "v4". Another processing will even- 
tually be performed using the final SI 8. 18.0 cali- 
bration files. New public versions of CASSIS will 
be posted as minor modifications in the algorithm 
warrant a new release, and the online database will 
be updated accordingly. 

3.1. Bad pixels 

The BMASK image contains the flag value of 
each pixel for a given exposure. CASSIS considers 
that a pixel is bad if the mask value is greater than 
256. Pixels with such flag values have the follow- 
ing issues identified by the SSC BCD pipeline: the 
flat field could not be applied, the stray-light re- 
moval or crosstalk correction could not be applied, 
pixel is saturated beyond correctable non-linearity 
in sample(s) along ramp, data missing in downlink 
in sample(s) along ramp, only one or no usable 
plane, pixel identified as permanently bad. We 
refer to the Spitzer/IRS documentation for more 
details on the flags. Although pixels with BMASK 
values above 256 could possibly be used, it is im- 
portant that the pipeline remains conservative in 
order to minimize potential artifacts. 

Another group of bad pixels, referred to as 
"rogue" pixels, are pixels that misbehave over long 
periods of time, randomly changing sensitivity on 
short timescales. They cannot be properly cal- 
ibrated. The SSC has released a series of masks 
identifying the long-term rogues for each IRS cam- 
paign. All the bad pixels, including rogue pixels, 
are cleaned using the IRSCLEAN tooll. The bad- 
fix method is used with the BMASK file along with 
a conservative "super-rogue" mask that combines 



The IRSCLEAN package can be found at 
|http://ssc.spitzer. caltech.edu/dataanalysistools/tools [ 
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Fig. 1. — Basic steps of the extraction pipeline. Only the steps leading to the optimal extraction are shown. 
Tapered column extraction are performed on the background-subtracted images, and the spectra are then 
defringed and nod-combined. The boxes with a dark background correspond to the reference branch which 
produces the best results in the majority of cases. The naming convention for each available product is 
shown in the bottom line in each box. 
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Fig. 2. — Basic steps of the extraction pipeline starting from the flux calibrated spectra. See Figure[T] for 
the figure description. 
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rogue pixels from the relevant observation cam- 
paign and earlier campaigns. The uncertainty file 
is cleaned using the same mask as for the data im- 
age. We refer to the IRSCLEAN documentation 
for more details on the cleaning process. 

Note that when several bad pixels are contigu- 
ous, the pixel replacement algorithm badfix will 
not fix the pixel(s) in the middle of the cluster. 
Further pixel verifications are performed in the 
next CASSIS pipeline steps to attempt a correc- 
tion. 

3.2. Co-addition of individual exposures 

Before co-adding the individual exposure im- 
ages (DCEs), CASSIS checks the dispersion of the 
reconstructed coordinates over the observation du- 
ratiorll. If the dispersion is significant, co- adding 
the images would result in blurring the source spa- 
tial profile so that the (optimal) extraction could 
become unreliable. For this reason, CASSIS co- 
adds the exposures only when the coordinate dis- 
persion is lower than a certain threshold. 

For a given observation, several cases can be 
distinguished for the image co-addition depending 
on the number of exposures available. 

1 exposure. Single exposures are simply 
transferred to the next step. 

2 exposures. Rather than performing a sim- 
ple combination (such as an error-weighted av- 
erage), the pipeline takes advantage of having 2 
exposures to produce a better result by fiagging 
bad pixels that were not cleaned or could not be 
cleaned fSection lXTj) . First, pixels are compared 
using their BMASK value and if one pixel has a 
higher value than the other, it is ignored. Further 
bad pixel flagging is achieved by analyzing each 
column of the image separately. A column can be 
seen almost as a spectrum, since the wavelength 
axis is almost parallel to the detector y-axis. Out- 
liers are then identified using their deviation from 
the local flux variance. When the 2 pixels are 
equally bad, the average is taken. In other cases 
when one pixel dominates the difference to the lo- 
cal variance, only the pixel from the other image is 



This information is provided by the header l^eywords 
SIGRA and SIGDEC populated by the SSC pipeline 
(Sect.©. 



used. All the other pixels that were not flagged as 
outliers have error- weighted average fluxes. The 
final uncertainty on the combined pixel is the sum 
of the average of the individual errors and of the 
flux difference between the 2 pixels to account for 
the fact that root-mean square (RMS) errors in 
the input images might be underestimated in some 
cases. 

More than 2 exposures. For each pixel, 
the median over all the images is averaged with 
the error-weighted-average. This is done to min- 
imize systematic errors due to unreliable uncer- 
tainty values. Uncertainties are combined accord- 
ingly, i.e., taking the average between the error on 
the median flux (median absolute deviation) and 
the error on the average (see also Sect. [7]). 

In a small fraction of the observations (« 
0.35%) there is a significant dispersion of the 
pointing coordinates over the exposures (0.05 pix- 
els, i.e., w 0.1" in SL and « 0.25" in LL) that 
requires the pipeline to performs steps in a differ- 
ent order; 

• The low-level rogue pixels are removed using 
the background images fSection [3.3|) . 

• The individual exposures are extracted sep- 
arately (Section[S]) . 

• The spectra corresponding to the individ- 
ual exposures are combined, resulting in one 
spectrum per module, order, and nod posi- 
tion. The two nod spectra are then com- 
bined fSection l^X^ . 

At this point of the pipeline, there is one image 
per module, order and nod position. Note that for 
staring cluster observations, there can be several 
cycles with several exposures for a given position. 
In this case, the spectra corresponding to the var- 
ious cycles are first combined. The online docu- 
mentation provides specific help for staring cluster 
observations. 

3.3. Removal of the background emission 
and of low-level rogue pixels 

Although images were cleaned fSection [3?T|) . 
low-level rogue pixels could remain. It is possible 
to remove their contribution to the intrinsic source 
emission by subtracting a background image or a 
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set of images from the same module. This step 
also allows subtraction of the large-scale back- 
ground emission (mostly dominated by zodiacal 
dust emission). Two methods are used: 

Subtraction by nod. The (single) image cor- 
responding to the other nod position is used for 
the difference. The SMART- AdOpt algorithm is 
used to test the presence of a point-like source or 
partially extended source at the location of the 
current nod extraction. If the test is positive, 
another background subtraction has to be used. 
Note that very extended emission does not mod- 
ify significantly the source's profile, and it can be 
removed during extraction (Section[3]). 

Subtraction by order. The two nod im- 
ages corresponding to the other order are (error- 
weighted) averaged. Contamination at the loca- 
tion of the current nod extraction is also tested. 

When the image difference is performed, uncer- 
tainties are combined quadratically (see Sect. [7]). 
In some cases, there are contaminating sources in 
both the "by nod" and "by order" images, so that 
a third method is used, referred to as "in situ", 
which simply removes the extended emission dur- 
ing optimal extraction, i.e., without removing the 
low-level rogue pixels (Section[5|) . 

The final spectrum provided to the user al- 
ways corresponds to the background subtraction 
method leading to the best S/N (Sect.|8]). In 
most case, the best method is the subtraction 
"by nod" since it involves a background that is 
observed shortly before/after the science target 
(Sect. refsec:general). Note that the other version 
of the spectrum corresponding to the other back- 
ground subtraction method is available as an op- 
tional product. 

4. Pre-extraction diagnostics 
4.1. Spatial extent 

All the images (exposure-combined, background- 
subtracted) go through the SMART- AdOpt pro- 
gram, which first estimates the source spatial ex- 
tent. The source extent diagnostic is essential to 
determine which extraction method is suitable. 
For point-like sources, the preferred extraction is 
the optimal extraction which produces the best 



S/N (Sect. 15.1) ). For partially-extended sources, a 
"tapered column" extraction is preferred since it 
recovers most of the source's flux. This extrac- 
tion includes all flux within the diffraction-limited 
spatial extent of the spectrum, without weighting 
by the PSF. For very extended sources, a full slit 
extraction is appropriate. 

In order to estimate the spatial extent, AdOpt 
derives the ratio between the full width at half 
maximum (FWHM) of the source spatial proflle 
and of the PSF profile. Both profiles are col- 
lapsed on the 20 detector rows corresponding to 
the shortest wavelengths, i.e., where the PSF is 
the narrowest. The resulting ratio is directly re- 
lated to the intrinsic spatial extent of the source: 

S = Py/x^ - 1, (1) 

where S is the intrinsic FWHM of the source, P 
is the FWHM of the PSF, and x is the ratio be- 
tween the FWHM of the observed profile and the 
FWHM of the PSF. Note that this method as- 
sumes that the source profile can be reproduced 
by a broadened PSF so it is accurate only for par- 
tially extended sources. 

The global extent of a source is calculated us- 
ing the average of the S determination for each 
module and order, assuming the extent is not a 
function of wavelength. Weights are applied to 
each value based on the detection level (Section[5|). 
FigureO gives an example of the output plot pro- 
duced to compute the spatial extent. 

Determinations of S can be different for each 
module because the intrinsic source extent might 
vary with wavelength. Further iterations of the 
CASSIS processing will include the calculation of 
the spatial extent for each wavelength element 

4.2. Multiple sources 

Before the extraction of the intended source is 
performed, the presence of another source in the 
slit is checked. Depending on the relative bright- 
ness between the intended source and the contam- 
inated source, it is possible that the local back- 
ground emission is not well determined (Sect.[5|). 
If there is a positive detection of a contaminat- 
ing source within the slit, the extraction continues 
but a flag is carried through for the corresponding 
module/order spectrum and appears in the result 
page of the online interface (Sect. [9]). In such cases 
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Fig. 3. — The intrinsic source extent S is calcu- 
lated for each module/order spectrum. Two mea- 
surements of S are available for each module and 
order, one for each nod position. The global extent 
is calculated using the average value of S weighted 
by the detection level in each spectrum. 

the users can check the diagnostics plots provided 
by the interface to judge the quality of the extrac- 
tion. The spectrum of the serendipitous source is 
not extracted. 

5. Spectral extraction 

Regardless of the source extent, both optimal 
extraction and tapered column extraction are per- 
formed. Hence, for a given source, there are 12 dif- 
ferent spectra (2 nod positions, 2 extraction tech- 
niques, 3 background-subtraction methods). De- 
pending on the spatial extent (Sect. l4TT|) . and on 
the best background available fSect. 13.31 ). the user 
will be presented with the best extraction choice 
and the other spectra will remain available as op- 
tional products. 

5.1. Optimal extraction 

Optimal extraction uses the PSF profile to 
weight the pixels from the spatial profile, while 
tapered column extraction integrates the flux in 
a spectral window that expands with wavelength. 
Optimal extraction is performed with the AdOpt 
algorithm which makes use of a super-sampled 
PSF (Lebouteiller et al. 2010). A super-sampled 
PSF is critical since it can be shifted and resam- 
pled anywhere along the slit, making the optimal 
extraction valid for any source's position. For this 



reason the AdOpt optimal extraction is an ideal 
choice for extracting the full set of IRS observa- 
tions. 

Since the CASSIS atlas is initially meant to pro- 
vide the spectra of targeted sources, the source 
finder is limited to positions around a given nod, 
with a range of ±2 px (see the pixel size in Tabled]) 
to account for slight mispointings. The range 
around the nod position ensures that the right 
source is extracted even when there is a brighter 
source in the slit. The source position is then fine- 
tuned to an accuracy of better than a tenth of a 
pixel. Examples of spectra are shown in FigurelD 

Depending on the geometry, extended emission 
might still be present despite the background im- 
age subtraction (Sect. 13. 3t . This is because the 
background image does not correspond exactly to 
the background at the source's position. For this 
reason, a 0-degree polynomial background is cal- 
culated on-the-fly by AdOpt. It is important that 
other bright sources are not contaminating the ob- 
servation so that the polynomial background is 
well estimated (Section|4|). 

In order to calculate the flux at a given wave- 
length, SMART- AdOpt performs a multiple linear 
regression (Lebouteiller et al. 2010). A quanti- 
tative detection level is calculated by the AdOpt 
algorithm. The detection level represents the per- 
centage improvement between the initial image 
and the residual image in which the source has 
been subtracted. For simplicity, the detection 
level is encoded into an integer value illustrating 
the quality of the detection: for sources not de- 
tected, 1 for tentatively detected, 2 for barely de- 
tected, 3 for detected, 4 for well detected. Note 
that the detection level is accurate only for point- 
like sources for which AdOpt is currently designed. 
For extended sources, a significant residual emis- 
sion can remain after the point-like source extrac- 
tion, effectively decreasing the detection level pa- 
rameter. 

5.1.1. Wavelength grid 

As explained in Lebouteiller et al. (2010), 
AdOpt extracts detector rows instead of pseudo- 
rectangles, the latter corresponding to a zone in 
the image where pixels have the same wavelength. 
As a result, the output wavelength grid of the 
extracted spectrum depends on the exact source 
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A0RKEY=28 1 1 1 872 (v4) (nods, no defringing, full sampling) 

{NED= NGC 6240:[GVA2004] NT (multiple)} [SIMBAD= [PMS2007] 13 (multiple)} 
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Fig. 4. — Example of a CASSIS spectrum for a given source. Only the optimal extraction of the spectrum is 
shown since the source is point-like. The tapered column extraction always remain available as an optional 
product. Top — The spectrum at each nod position is plotted and the ratio between both spectra is given 
in the bottom panel. This is an optional plot that the user can choose to inspect for diagnostic purposes. 
Bottom — The two nod spectra have been combined, providing the default product. Vertical dotted lines 
indicate the position of the brightest atomic and molecular lines. This is the main plot output by the web 
interface. 
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position. Optionally, the wavelength grid can be 
interpolated afterwards on a common reference 
grid (the SSC "wavesamp" calibration file). 

The choice of the wavelength grid depends 
whether spectral resolution or S/N should be priv- 
ileged. By choosing the observed wavelength grid, 
the spectral resolution is slightly better than what 
can be achieved with the optimal extraction of 
the SSC SPICE software at the expense of S/N. 
In contrast, the S/N when interpolating the wave- 
length grid is similar to the optimal extraction in 
SPICE at the expense of spectral resolution. It 
must be kept in mind however that the wavelength 
grid interpolation results in a smoother spectrum 
but it can make two bad pixels out of one (which 
is also an undesired effect of pseudo-rectangle ex- 
tractions). 

Since the choice of the wavelength grid is not 
only a technical but also a scientific decision, the 
CASSIS interface always proposes both options to 
the user. By default however, CASSIS uses the in- 
terpolated spectra in order to provide wavelength 
grids compatible with the tapered column extrac- 
tions. 

5.1.2. Nod spectra combination 

Although the spectra of each nod are available 
as optional products, CASSIS produces the nod- 
combined spectrum, which is the default output. 
The combination process depends on whether the 
spectra were interpolated or not (Section lS.l.ip . 
We discuss both methods separately. The error 
treatment is described in Section[71 

Reference "wavesamp" wavelength grid 
(standard product provided as the default). 
The two nod spectra are calculated on the same 
reference wavelength grid. The co-added spec- 
trum is the error-weighted average of the two nod 
spectra. The error function (difference between 
the spectra) is then calculated, smoothed using 
a multi-resolution algorithm, and used to identify 
outliers in the individual nod spectra. Pixels are 
corrected accordingly to their relative discrepancy, 
i.e., if the pixels in the two nod spectra are outliers, 
the error- weighted average is used, but if only one 
pixel is an outlier, the other nod spectrum is used. 

Observed wavelength grid (optional prod- 
uct). The two nod spectra are first interleaved 
and aligned (see Figure 4 in Lebouteiller et al. 



2010). Alignment is performed by calculating the 
smoothed error function (difference between the 
spectra). The result is equally split to each nod 
to align the spectra. Outliers are identified es- 
sentially the same way as for the interpolated 
nod spectra. The only difference is that a pixel 
identified as an outlier is flagged and given a not 
a number (NaN) value. The final spectrum has 
twice as many points compared to the individual 
input nod spectra. It is referred to as the "fully 
sampled" spectrum. 

Another optional version of the final spec- 
trum is calculated by interpolating the interleaved 
spectrum on the reference wavelength grid. Note 
that this is different from the co-addition of in- 
terpolated nod spectra. The interpolation of the 
interleaved spectrum includes more points and is 
usually more accurate. 

5.1.3. Flux calibration 

Flux density is converted from e-/sec to Jy us- 
ing the default option in SMART- AdOpt, i.e., the 
use of a relative spectral response function (RSRF) 
derived from the comparison between observations 
of calibration stars and their theoretical models. 
The calibration star HR6348 was used to pro- 
duced the CASSIS flux calibration for all mod- 
ules except LL after campaign 45 (the LL detector 
bias changed after this campaign). To calibrate 
LL from campaign 45 and on, RSRFs were con- 
structed using the stars HR 6348 and HD 173511. 
A future paper (G.C. Sloan et al., in preparation) 
will explain the construction of the templates for 
these two stars and their photometric calibration. 

For point-like sources, the CASSIS inter- 
face displays by default the spectrum extracted 
with optimal extractior{^. For partially-extended 
sources, the CASSIS interface displays the ta- 
pered column extraction, which provides a better 
flux calibration although there is still no satisfac- 
tory means of accurately calibrating such sources 
fSection [5?2|) . Note that an upgrade of SMART- 
AdOpt is underway to perform optimal extrac- 
tions of slightly extended sources. The CASSIS 



Using optimal extraction on a partially extended source 
would result in a slightly underestimated flux density and 
would slightly modify the spectrum slope, an effect which 
also exists, though to a lesser extent, for a regular tapered 
column extraction. 
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database will eventually include these spectra with 
the corresponding calibration. 

In about 3% of the observations the source is 
not well centered in the slit in the dispersion di- 
rection, requiring a special flux calibration. In 
such cases, a fraction of the PSF lies outside the 
slit so that the regular optimal extraction fails 
to fit the proper profile, providing an underesti- 
mated fiux. Although the manual extraction of 
AdOpt can solve this problem by modifying on- 
the-fly the PSF profile as a function of the shift 
in the dispersion direction, this is not possible 
automatically. Instead, CASSIS checks the PT- 
GDIFFY header keywortS which gives the point- 
ing error in the dispersion direction between the 
coordinates requested by the user and the field of 
view coordinates (effective coordinates at the nod 
position in the current module and order). This 
is by no means a definite proof of a genuine off- 
set as the MIR centroid might not coincide with 
the requested coordinates. In any case, the data is 
flagged if PTGDIFFY is larger than a given frac- 
tion of the slit height and the user is then advised 
to perform a manual extraction. 

5.2. Tapered column extraction 

For partially extended sources, the default out- 
put spectrum is extracted with a tapered column. 
The tapered column extraction uses an extraction 
aperture whose width varies with wavelength be- 
cause of the varying diffraction limit with wave- 
length. Depending on the spatial extent derived 
in (Sect. HI, the extraction aperture width is scaled 
up to account for all of the source flux. 

The flux calibration for tapered column extrac- 
tion uses the default SSC flux conversion tables. 
Note that while the flux might be closer to the 
truth when using tapered column extraction on 
a partially extended source as compared to using 
optimal extraction, there is currently no satisfac- 
tory calibration for sources other than point-like 
sources and very extended sources (significantly 
larger than the slit). One undesirable effect is 
that, when using a point-like source flux calibra- 
tion (whether in optimal extraction or tapered col- 
umn extraction) on a partially extended source, 
the LL spectrum contain more flux than the SL 



spectrum spectra because the LL extraction aper- 
ture is larger than the SL aperture. A more ad- 
vanced version of the flux calibration is currently 
under investigation and will be included in future 
CASSIS processings. It uses a theoretical PSF and 
computes the amount of light lost outside the slit 
as a function of a broadening parameter. 

6. Defringing 

Fringes are common in infrared detectors, they 
produce sinusoidal variations in the spectrum 
(e.g., Kester et al. 2003). For Spitzer/IRS, spectra 
from the LLl module are the most affected. The 
IRSFRINGElil] algorithm is used by CASSIS to 
remove the fringes in that module only. We lim- 
ited the defringing solution to two sine functions 
to prevent an overcorrection of any real structure 
in the spectrum. Furthermore, since defringing 
can produce undesired artifacts for sources with 
a low S/N, only sources with LLl spectra with 
S/N^5 are considered. 

Note that since defringing is a complex process 
that uses several hypotheses and parameters, we 
also provide the uncorrected spectra available as 
an optional product. 

7. Flux Uncertainties and Calibration 

For all the steps until the spectra combination 
(Sect. [5T!2|) . the statistical (RMS) errors within 
an individual spectrum are fully propagated using 
the standard equation: 



E^=0 ^'i'^f 



N 



(2) 



where w is the weight and a is the uncertainty 
for the image (or spectra) number i. The weight- 
ing factor, although arbitrary, depends on the er- 
ror distribution. The error distribution is mostly 
poissonian (photon count noise) for sources with 
fluxes signiflcantly above the detector readout 
noise, while it becomes "normal" for very faint 
sources. The pipeline processing steps then nat- 
urally tend to "normalize" errors. We chose to 
consider weights as l/cr^ throughout the pipeline. 



^"The value of this keyword is calculated by the SSC BCD 
pipeline (Sect.[2]|. 



"IRSFRINGE 



be 



found 



at 



_http: / /ssc.spitzer.caltech.edu/dataanalysistools/tools. 
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The combination of the spectra introduces a 
systematic error illustrated by the flux density dif- 
ferences between the two nod spectra. This differ- 
ence is usually due to the presence of other sources 
in the slit affecting the background subtraction, or 
to pointing errors resulting in a slight shift in the 
dispersion direction. The two errors (RMS and 
systematic) are given as two different fields in the 
output files. 

Another source of systematic error is the fiux 
calibration itself (G.C. Sloan et al. in prepara- 
tion). The corresponding uncertainty depends on 
the stellar template(s) used (i.e. on the uncer- 
tainties inherent to the stellar models themselves) 
and on the dispersion in the various spectra of the 
calibration stars used to build the RSRF. With 
the current flux calibration in CASSIS, we esti- 
mate this systematic error to be less than « 2%. 
Wavelength ranges near the edges of detectors are 
not so well calibrated because of instrumental ef- 
fects such as light leaks. Such uncertainties are not 
yet estimated for the first public CASSIS version 
("version 4"). 

8. Products and diagnostics 

The main and final product is the calibrated 
spectrum resulting from merging the spectra of 
the various modules, nods, and spectral orders. 
The default spectrum presented to the user de- 
pends whether the source is point-like (optimal 
extraction was selected) or partially-extended (ta- 
pered column extraction was selected). The de- 
fault spectrum includes defringing (Sect. [6]), and 
the wavelength grid is the reference "wavesamp" 
(Sect. lS.lTTj) . The background subtraction method 
used for a given module and order (Sect. [231) is the 
background providing the best detection level ac- 
cording to the AdOpt algorithm. In less than 2% 
of the observations, the presence of a contaminat- 
ing source in the background image forces the de- 
fault spectrum to be presented without any back- 
ground subtraction. 

The other (non default) versions of the spec- 
tra are also available as optional products (op- 
timal or tapered extraction, specific background 
subtraction method, specific wavelength grid, no 
defringing, etc.). These versions are intended for 
advanced users who wish to inspect in detail the 
spectra and to compare the various extraction 



techniques or even the various wavelength scales 
available. The online documentation and inter- 
face advises users to check systematically the op- 
tional products especially in the first year of CAS- 
SIS since improvements might be needed early on. 
The various products and the naming convention 
are the ones shown in Figures[l] and [21 

The most important parameters defining the 
output diagnostics for a given AORkey are the 
coordinates, reconstructed from the source's po- 
sition inferred for each spectral extraction with 
AdOpt; the detection level, corresponding to the 
maximum detection level in the spectra of each 
nod, module, and order; and the spatial extent, 
described in Sect. 14. II 

In the final step of the pipeline, the extracted 
source is resolved using the Nasa/IPAC Extra- 
galactic Database (NElQ and the SIMBAD As- 
tronomical Databasj^. In cases where multi- 
ple sources fall within the PSF beam, the clos- 
est matching source is chosen and a special flag 
is given. The spectroscopic redshift is inferred 
using NED and the coordinates of the resolved 
source. No redshift is returned if multiple sources 
are found within the PSF beam. The source iden- 
tification and redshift are provided to the user 
and are part of the CASSIS atlas. Photometric 
and spectroscopic redshifts from VIZIElS cata- 
logs will be included in future CASSIS versions. 

9. Online interface 

The online CASSIS interface offers basic func- 
tionalities such as querying by AORkey, program 
ID, or coordinates. Eventually, additional param- 
eters will be added to allow querying in various 
other ways, such as by source spatial extent or 
spectral properties, or the redshift (Sect. [51). 

The results page includes observational param- 
eters, extraction diagnostics, and the links to 
download spectra. The interface also enables view- 
ing and analyzing the spectra with Virtual Obser- 
vatory (VO) tools. We currently interoperate with 
the International Virtual Observatory Alliance at 
two levels: we permit downloading of spectra in 



http://nedwww.ipac.caltech.edu/ 
http: / /simbad. u-strasbg.fr/simbad/ 
http: / /vizier. u-strasbg.fr/viz-bin/VizieR 
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the IVOA Spectral Data ModeS so that they 
may be processed with VO spectral analysis tools 
such as VOSpec and SPLAT. We also allow run- 
ning the VO-enabled SPLAT applet directly on a 
given dataset. Eventually, we plan to support Sim- 
ple Spectral Access Protocol (SSAF0) discovery 
through Virtual observatory search tools. 

10. CASSIS Source Statistics 

The CASSIS atlas currently includes 13 264 dil- 
ferent observations, i.e., most of the low-resolution 
observations ever performed in the staring mode 
(single and cluster mode) by the IRS instrument. 
Figure[S] shows the distribution of observations in 
the sky. Among the CASSIS sources, about 90% 
are at least barely detected. According to the 
detection parameter from AdOpt, we find that 
« 73% of the objects are well detected (level 4, 
see SectEj]), w 5% are detected (level 3), « 7% 
are faint (level 2), w 4% are barely detected (level 
1), and « 10% are not detected (level 0). 

Figure|6^ shows the histogram of the spatial ex- 
tent distribution. Most sources (70%) have an in- 
ferred intrinsic extent below 2" and can be consid- 
ered point-like (see Sect. 14. ip . FigurelHb shows the 
histogram of the pointing offsets, i.e., the source 
position in the cross-dispersion direction of the 
slit. Most sources are found within a fraction of a 
pixel from the exact nod position. 

For about 2118 distinct sources, a redshift 
could be associated using NED and the spatial 
coordinates of the extracted source (see Sect. [8]). 
The redshift distribution of the detected sources 
is shown in Figure[7l We selected the sources 
with significant PAH emission from this sample 
to build the stack of Figure[5] Spectra are plot- 
ted in the rest-frame. The good alignment of the 
known spectral features indicates that the redshift 
association provides reliable values. Note that the 
sources in this example span a wide variety of ob- 
jects, the only common property being a spectro- 
scopic redshift determination. 

11. "Build your own" samples 

A specific aspect of CASSIS is that the fiux den- 
sity at every wavelength element of the reference 




Fig. 6. — Top ~ Spatial extent of sources with a 
level detection above 2, i.e., sources at least barely 
detected (see Sect. lS.ll for the determination of the 
detection level). The vertical dotted line corre- 
sponds to an extent of 2" below which sources 
can be safely considered as point-like. Bottom — 
Pointing offset distribution (in the cross-dispersion 
direction) of sources with a level detection above 
2, i.e., sources at least barely detected. The 2 ver- 
tical dotted lines indicate a fraction of 0.2 pixels 
in SL (0.36") and in LL (1.2"). 



^^http:/ /www. ivoa.net/Documents/SpcctrumDM/index. html 
■'^http: / /www. ivoa.nct/Documents/latcst/SSA. html 
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Fig. 5. — IRS observations projected on the sky. The top figure shows the mapping observations and the 
moving objects which are not included in the atlas. The middle figure shows the observations with spectra 
extracted in CASSIS, corresponding to staring observations (single or cluster). The bottom figure shows the 
subset of sources from CASSIS with associated redshifts. 
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Fig. 8. — Rest-frame spectra of extragalactic objects are plotted as rows, the light scale scaling with the 
flux density. Only sources with significant PAH emission are included in this sub-sample. The redshift was 
associated with NED from the coordinates of the extracted source. A few black curves can be seen that 
corresponds to noisy regions in each module. 
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(wavesamp) table is contained in the database. 
Therefore, spectra themselves can be used in a 
query. For instance, it is possible to query the 
database for spectra matching a given reference 
spectrum (model, user-defined, or existing spec- 
trum in the database). The match can be per- 
formed on the full wavelength scale or on any sub- 
region, for example around a given spectral feature 
(see example of Figure|9]). 

Furthermore, it is possible to perform oper- 
ations on the spectra to satisfy a given query. 
For instance, a few simple examples of semantic 
searches could be: 

- "select sources where /(14,0.2) > /(5,0.2)", for 
spectra with a positive slope from 5 /im to 14 /im. 
The continuum flux is calculated by integrating a 
0.2 /im window around a central wavelength, 

- "select sources where /(6.2, 0.2)-0.5x [/(5.8, 0.2)-^ 
/(6.4,0.2)] > 0.1", for spectra with significant 
emission in the 6.2 /im PAH band. 

- "select sources where SNR(10, 1) > 5", for spec- 
tra in which the S/N is greater than 5 between 9 
and 11 /im. 

- "select sources where (extent< 2 and glohaldetlvl— 
4 and z < O.lj", for observations of point- like 
sources with a detection level equal to 4 (well de- 
tected) and with a redshift below 0.1. 

More sophisticated measurements (e.g., involv- 
ing line fits) for the full CASSIS sample will be 
provided in a later step. The sample of extra- 
galactic sources with known redshifts is currently 
being investigated, including the 6.2 /tm PAH and 
the silicate strength values (H. Spoon et al., in 
preparation). The results will be eventually in- 
corporated in the CASSIS database and will be 
usable by queries. The mid-infrared classification 
from Spoon et al. (2007), using these 2 spectral 
parameters will be also included in the database. 

A few illustrations of extragalactic templates 
built using spectra in the CASSIS database are 
shown in Figures[TO] and[TlJ These spectra can be 
considered as representative of the objects show- 
ing a particular spectral feature. The template 
shows the median spectrum for the 100 sources 
with strongest features using a given criterion. 
In these examples, emission from lines and PAHs 
were simply measured by integrating the flux, no 
fit was performed. The silicate absorption opti- 
cal depth was measured according to Spoon et al. 



(2007), assuming the method usually applicable to 
PAH-dominated spectra. There is no selection ac- 
cording to object type, so the spectra include the 
contribution from starburst galaxies, active galac- 
tic nuclei, and all other sources in the CASSIS 
atlas. 

Figure[TO] illustrates preliminary investigations 
for correlating the presence of a given line or band 
with other spectral features. The PAH-dominated 
template agrees well with the starburst template 
spectra from Brandl et al. (2006) and Bernard- 
Salas et al. (2009). Deeper silicate absorption 
is observed in our template, which is likely due 
to the contribution from ultraluminous infrared 
galaxies in addition to starburst galaxies. An- 
other interesting result is the correlation between 
H2 emission and the [Si 11] emission at 34.82 /im 
(e.g., Roussel et al. 2007). The objects showing an 
intense [Fe ii] + [0 iv] line complex tend to show 
strong [Ne ill] emission at 15.56 /im. Finally, the 
objects with the deepest silicate absorption show 
few lines and weak PAH bands. 

Figure[Tl] splits the sample of PAH-dominated 
sources with the most extreme 6.2 /tm/ 11.3 /im 
PAH ratios. It can be seen that objects with a 
large 6.2/11.3 ratio show a steeper dust continuum 
while also having a larger [S iv]/[Ne 11] ratio, sug- 
gesting that the PAH band ratio correlates with 
the radiation field hardness (e.g., Lebouteiller et 
al. 2011; Madden et al. 2006; Wu et al. 2006; En- 
gelbracht et al. 2005). 

12. Summary 

We present the CASSIS atlas of Spitzer IRS 
low resolution spectra. All IRS observations in 
staring mode are included, totaling about 13 000 
spectra corresponding to 11000 distinct sources. 
The output spectra provide the best product for a 
publishable, measurable IRS spectrum. 

An online interface (http://cassis.astro.cornell.edu/atlasl 
accesses the spectra, provides detailed information 
on production of the spectra, and allows searches 
of the atlas based on various parameters. 

Two versions of the spectra are available, using 
the optimal extraction provided by AdOpt (suited 
for point-like sources) and the tapered column 
extraction (better suited for partially extended 
sources). Several important diagnostics are pro- 
vided, most notably a quantitative detection level 
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Fig. 9. — The red spectrum (corresponding to a star) was given as an input to find matching spectra in the 
database. The criterion in this case was to match the spectra over the whole wavelength range. The first 
hit (left) corresponds to the input object itself. The few first other hits are shown on the right. The match 
can be performed over any wavelength range (observed or rest-frame) . 
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Fig. 10. — Average spectra of cxtragalactic objects showing significant 6.2 /im PAH emission (a), 17.05 /Ltm 
H2 emission (b), [Fe ii]+[0 iv] (blended lines) emission (c), and silicate absorption optical depth (d). In 
each example, the 100 most extreme objects were selected. 
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Fig. 11. Average spectra of extragalactic objects showing the highest and lowest 6.2/um/11.3/im PAH 
band ratio. The 100 most extreme cases were selected in each case. 
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and the spatial source extent. 

Future versions of CASSIS will include point- 
like sources observed serendipitously (at other po- 
sitions in the slit or within the offset slit), sources 
observed with spectral mappings, and the IRS 
high-resolution observations. 

We wish to thank the people who contributed 
to the data reduction efforts over the IRS mis- 
sion. Former ISC members are especially ac- 
knowledged (in particular D. Devost, D. Levitan, 
D.Whelan, K. Uchida, J.D.Smith, E.Furlan, 
M. Devost, Y.Wu, L.Hao, B.Brandl, S.J.U. Hig- 
don, P. Hah) for their work on the SMART 
software and for the development of reduction 
techniques. Moreover, the Rochester group (in 
particular M. McCLure, C. Tayrien, I. Remming, 
D.Watson, and W.Forrest) played an important 
role in collaborating with the ISC to bring addi- 
tional and essential improvements to the data re- 
duction used in CASSIS. This research has made 
use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and 
Space Administration. This research has made 
use of the SIMBAD database, operated at CDS, 
Strasbourg, France. This research was conducted 
with support from the NASA Astrophysics and 
Data Analysis Program (Grant NNX10AD61G). 
Finally, we thank the referee for his/her useful 
comments. 
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Fig. 7. — Redshift distribution of sources with a 
level detection above 2, i.e., sources at least barely 
detected (see Sect. lS.ll for the determination of the 
detection level). 
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Table 1: Modules of the Spitzer/IRS. CASSIS currently includes only the data from the SL and LL modules. 



Module 


Order (s) 


A ( /im) 


Aperture size (") 


Pixel size 


(") 


SL 


1 


7.4 - 14.5 


3.7x57 


1.8 




SL 


2 


5.2 - 7.7 


3.6x57 


1.8 




SL 


3 


7.3 - 8.7 


3.6x57 


1.8 




LL 


1 


19.5 - 38.0 


10.7x168 


5.1 




LL 


2 


14.0 - 21.3 


10.5x168 


5.1 




LL 


3 


19.4 - 21.7 


10.5x168 


5.1 




SH 


11-20 


9.9 - 19.6 


4.7x11.3 


2.3 




LH 


11-20 


18.7 - 37.2 


11.1x22.3 


4.5 





Table 2: Number of observations and distinct sources in the CASSIS atlas per scientific category. 



Category 


Observations^ 


Detections'^ 


Distinct sources'^ 




1252 


504 




Intermediate-z galaxies 


847 


753 


753 


Nearby galaxies 


484 


289 


286 


Local Group galaxies 


561 


546 


542 


AGN / quasars /radio-galaxies 


1535 


1424 


1385 


Interacting/ mergers 


156 


147 


147 


ULIRGS/LIRGS 


631 


601 


601 


Starburst galaxies 


123 


120 


120 


Galaxy clusters 


57 


52 


51 


Cosmology 


3 


1 


1 


Extragalactic jets 


5 


5 


2 


Gamma-ray bursts 


5 


5 


5 


Compact objects 


57 


56 


53 


Galactic structures 


13 


12 


12 


Massive stars 


155 


118 


118 


Evolved stars 


1256 


1164 


1074 


Brown dwarfs 


287 


282 


274 


Stellar population 


218 


218 


164 


Star formation 


577 


548 


526 


Young stellar objects 


1670 


1556 


1492 


Circumstellar disks 


2464 


2417 


2379 


Extra-solar planets 


7 


7 


4 


ISM 


884 


768 


768 


HII regions 


57 


44 


44 


Total 


13264 


11637 


11304 



Note. — The following categories currently have no spectra in CASSIS: satellites, comets, planets, kuiper-belt objects, near- 
earth objects, asteroids, star clusters, zodia<;al dust, and dark matter. 
"The number of observations includes the number of AORkeys and the number of pointings for the cluster observations. 
''Detection level greater than 1, i.e., sources at least barely detected. The difference between the number of observations and the 

number of detections is often due to non-detections in background observations designed as separate AORkeys. 
'^Distinct sources are found by matching observations within 4" of each other. 
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Appendices 



A. Acronyms and abbreviations 

AdOpt - "Advanced Optimal extraction". Plugin program in the SMART environment notably enabling 
optimal extraction of Spitzer/IRS spectra by using a super-sampled point spread function. 
AORkey (or AOR) - "Astronomical Observation Request key". Unique identifier for observations per- 
formed by Spitzer. There can be several AORkeys per object. 

BCD - "Basic Calibrated Product". Calibrated IRS detector images (128 x 128 pixels) provided by the SSC 
pipeline. 

BMASK - "Bit-mask". Image plane (128 x 128 pixels) containing a possible error condition code for each 

pixel in a detector of the IRS. The bit codes can be found in this page: http://isc.astro.cornell.edu/SmartDoc/ErrorProcessing 

DCE - "Data Collection Event". Single exposure images also referred to as "ramps". The ramp duration is 

the time between the first and last non-destructive reads of the array. Several ramp times are available for 

each module. 

IRS - "InfraRed Spectrograph" . The spectrograph onboard the Spitzer Space Telescope. 

LL - "Long-Low" . One of the low spectral resolution modules of the IRS (Table[T]) . 

PSF - "Point Spread Function" . Instrument spatial response function to a point-like source. 

RSRF - "Relative Spectral Response Function". Function giving the ratio of the observed spectra of a 

calibration star in uncalibrated units to the theoretical stellar template. 

SMART - "Spectroscopic Modeling Analysis and Reduction Tool". Contributed software allowing the 
community to reduce and analyze data from the IRS. 

SPICE - "SPitzer IRS Custom Extractor" . Official software developed by the SSC to reduce IRS data. 
SSC - "Spitzer Science Center" . Organization at CalTech supporting the science for Spitzer. 
SL - "Short-Low". One of the low spectral resolution modules of the IRS (Table[T|). 
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